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^ ■ Abstract 

m ! 

, We continue our previous application of supersymmetric quantum mechanical 

' methods to eigenvalue problems in the context of some deformed canonical commuta- 

tion relations leading to nonzero minimal uncertainties in position and/or momentum. 
fH ! Here we determine for the first time the spectrum and the eigenvectors of a one- 

dimensional harmonic oscillator in the presence of a uniform electric field in terms 
of the deforming parameters a, (3. We establish that whenever there is a nonzero 
minimal uncertainty in momentum, i.e., for a 7^ 0, the correction to the harmonic 
oscillator eigenvalues due to the electric field is level dependent. In the opposite case, 
^ ' i.e., for a = 0, we recover the conventional quantum mechanical picture of an overall 

. energy-spectrum shift even when there is a nonzero minimum uncertainty in position, 

i.e., for /3 7^ 0. Then we consider the problem of a i^-dimensional harmonic oscillator 
in the case of isotropic nonzero minimal uncertainties in the position coordinates, 
depending on two parameters (5, (3' . We extend our methods to deal with the corre- 
sponding radial equation in the momentum representation and rederive in a simple 
way both the spectrum and the momentum radial wave functions previously found by 
solving the differential equation. This opens the way to solving new D-dimensional 
problems. 
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1 Introduction 



During recent years, there has been much interest in studying quantum mechanical problems 
under the assumption of a modified Heisenberg uncertainty relation leading to nonzero 
minimal uncertainties in position and/or momentum. Such works are motivated by several 
independent lines of investigations in string theory and quantum gravity, which suggest 
the existence of a finite lower bound to the possible resolution of length Axq (see, e.g., 
[1, 2, 3]). Furthermore, the absence of plane waves or momentum eigenvectors on generic 
curved spaces also hint at a finite lower bound to the possible resolution of momentum Apo 
(see, e.g., [4]). 

Such nonzero minimal uncertainties in position and momentum can be described in the 
framework of small quadratic corrections to the canonical commutation relations [5, 6, 7]. 
These corrections can also provide an effective description of non-pointlike particles, such 
as quasiparticles and various collective excitations in solids, or composite particles, such as 
nucleons and nuclei [6]. 

The resolution of quantum mechanical problems with such deformed canonical commu- 
tation relations has been mostly restricted to the case where there is only a nonzero minimal 
uncertainty in position. Then one may indeed consider a deformed Schrodingcr equation 
in momentum representation and solve it using the technique of differential equations. An 
exact solution to the one-dimensional harmonic oscillator problem has been obtained in this 
way [8]. This approach has been extended to D dimensions [9] and some ladder operators 
have been constructed [10]. Some perturbative or partial results have also been obtained 
for the hydrogen atom [11, 12]. 

The case where there are nonzero minimal uncertainties in both position and momentum 
is much more involved because there is neither position nor momentum representation, so 
that one has to resort to a generalized Fock space representation or, equivalently, to the 
corresponding Bargmann representation [5, 13]. In this context, we have recently solved [14] 
the eigenvalue problem for the one-dimensional harmonic oscillator in a purely algebraic way 
by availing ourselves of an extension of supersymmetric quantum mechanical (SUSYQM) 
and shape-invariance techniques. The SUSYQM formalism [15, 16] supplemented with 
shape invariance under parameter translation [17, 18] is known to be a reformulation of 
the factorization method developed by Schrodinger [19] and by Infeld and Hull [20] (for 
a comparison between these two methods and corresponding references see [16, 21]). The 
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procedure that we have used rehes on a generahzed type of shape invariance, namely that 

connected with parameter scaling [22, 23, 24, 25]. 

The purpose of the present paper is to further illustrate the power of SUSYQM tech- 
niques in solving eigenvalue problems corresponding to deformed canonical commutation 
relations. 

First, we will consider the case of a one-dimensional harmonic oscillator with nonzero 
minimal uncertainties in position and momentum in the presence of a uniform electric field. 
Contrary to the solution of the corresponding problem in conventional quantum mechanics, 
which can be obtained from the solution without electric field by a simple coordinate shift, 
that of the present problem is more involved and will reveal some new features. 

Next, we will show how to extend our method to higher-dimensional problems by pro- 
viding an alternative and simple solution to the eigenvalue problem for the D-dimensional 
harmonic oscillator with isotropic nonzero minimal uncertainties in the position coordinates, 
which was previously dealt with by the differential equation technique [9]. 

Our paper is organized as follows. The one-dimensional harmonic oscillator in a uniform 
electric field is considered in section 2. Section 3 deals with the D-dimensional harmonic 
oscillator. Finally, section 4 contains the conclusion. 

2 One-dimensional harmonic oscillator in a uniform 
electric field 

Let us consider the case of a particle of mass m and charge ^ in a harmonic potential ^mcu^x^ 
and a uniform electric field £ parallel to the x-axis. It is described by the Hamiltonian 

1 - 

H + -mc;V - q£x (2.1) 

2m 2 ^ ^ ^ 

and the corresponding eigenvalue problem reads 

= £;„|Vn) n = 0,l,2,.... (2.2) 

Here x and p are assumed to satisfy the deformed canonical commutation relation 

[x,p] = ih{l + ax^ + Pp^) (2.3) 

where a > 0, (3 > 0, and a(3 < ft^^, so that the minimal uncertainties in position and mo- 
mentum are given by Axq — h^J(3/{l — h'^a(3) and Apo = h^a/ (1 — h^afi), respectively [5] . 
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In terms of dimensionless operators, X = x/a, P = pa/h, h = H/ {hw) and dimensionless 
parameters a = aa? and [3 = Ph^/a^, E = q£a/ (hu), where a = ^Jh/ (ma;), equations (2.1) - 
(2.3) can be rewritten as 

h^hp^ + X^)-SX (2.4) 
^iV'n) = e„|'0n) en = Ej{huj) n = 0,1,2,... (2.5) 
[X, P] = i(l + aX"" + /3P2) (2.6) 

with a > 0, /3 > 0, and q;/3 < 1. 
2.1 Energy spectrum 

At this stage, it is worth noting that going to transformed operators X' = X — £ and 
P' = P, as in conventional quantum mechanics, would convert h into a shifted harmonic 
oscillator Hamiltonian 

h^hp'^ + X'^)-ls^ (2.7) 
but, at the same time, change the commutation relation (2.6) into 

[X', P'] = i(l + aX'^ + (5P'^ + 2a£X' + aS''). (2.8) 

Hence, except for a = 0, the energy spectrum of h cannot be deduced from that for £ = 
obtained in our previous paper [14] (henceforth referred to as I and whose equations will 
be subsequently quoted by their number preceded by I). 
Let us instead try to factorize h as 

h^B+{g,s,r)B-{g,s,r) + e^ (2.9) 

where 

B^{g, s, r) = i=(sX T igP + r) (2.10) 

and eo is the factorization energy. In (2.10), g, s and r are assumed to be three constants 
depending on the parameters a, P, £ of the problem. The first two, g and s, are chosen 
positive and going to 1 in the hmit a, (3 0, while r is taken as going to —£ in the 
same hmit. With this choice, the operators B^{g, s,r) and B~{g, s,r) will be counterparts 
of the creation and annihilation operators for the shifted harmonic oscillator occurring in 
conventional quantum mechanics. 
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On inserting (2.10) in (2.9) and comparing the resulting expression with (2.4), we get 
four conditions, of which the first two are independent of 8 and coincide with equations 
(12.7) and (12.8). They fix the values of g and s, 

g^sk s=^=L= k=\{(5-a) + ^l + \{(5-ay (2.11) 
Vl — 

in terms of a and 13. The remaining two conditions read 

rs = -£ (2.12) 

= ligs-r'). (2.13) 

Equation (2.12) provides the value of r, 

r = -£Vl - ak (2.14) 

while equation (2.13) leads to the value of the factorization energy in terms of the three 
parameters a, (3, 8. 

Having proved that the Hamiltonian h is factorizable, let us now consider a hierarchy 
of Hamiltonians 

i 

hi = B+{gi,Si,ri)B-{gi,Si,ri) + ^ej i = 0,1,2,... (2.15) 

j=0 

whose first member Hq coincides with h. Here g^, Si, e^, are some parameters, the first 
three being assumed positive, and go = g, sq = s, tq = r. 

Proceeding as in I, let us impose the shape invariance condition 

B^{gi, Si, ri)B^{gi, Si, n) = B^{gi+i, Si+i, ri+i)B~{gi+i, s^+i, r^+i) + e^+i (2.16) 

where i = 0, 1, 2, It is equivalent to a set of four relations, of which the first two are 

again independent of £ and coincide with equations (12.15) and (12.16), respectively. In I, 
the latter have been solved by introducing some new combinations of parameters 

Ui = gi + ISi gi- -fSi 7 = J- (2.17) 

V a 

V ■ 

di = UiVi ti = — (2.18) 

Ui 

thereby leading to 

di^d ti^ q-H or Ui = q'/'^u Vi = q'^^'^v (2.19) 



where 

d = uv t = - q = ^^=. 2.20 

u 1 — y/ap 

The remaining two relations read 

Tj+iSi+i = riSi (2.21) 

On taking equations (2.12), (2.17), (2.19), and (2.20) into account, we obtain from (2.21) 
the solution for rj, 

--^'-'l-?)"- ''''' 

Furthermore, equation (2.22) leads to the eigenvalues e„ of h, 

n n—1 

en{q, t,S) = ^ei = Y^ QiSi + -QnSn - -T^ = e„(g, t, 0) + Ae„(g, t, 8) (2.24) 

1=0 1=0 ^ ^ 

where 



(2.25) 

are the eigenvalues in the absence of electric field, given in (12.30), and 



i^en{q. t, £) = ^g-" I 1 - - 1 (2-26) 



are the corrections due to the electric field. The latter can be rewritten as 

Ae„(g, z) = -\k'z\1 - tfq-^ (l - ^ (2.27) 
in terms of a new parameter proportional to 



z = A'yJ ^ ^ , ,5 — r (2.28) 

which will subsequently prove convenient (see equations (2.47) and (2.57)). 

In contrast with conventional quantum mechanics, the additional negative contributions 
Aen{q,t, z) are n dependent for generic a and f3 values. Since q > 1, their absolute value 
decreases from ^K^z^ to when n goes from to 00. As it can be inferred from the 
definition (2.4) of h, where the electric field breaks the symmetry under exchange of X and 



P, /S.en{q-it,z) is not invariant under exchange of a and (3 as it is the case for e„(g,t,0). 
This will lead to different hmiting cases in the next subsection. 

The corrections to the excitation energies in the absence of electric field 

en{q, t, 0) - eo(q, t, 0) = ^K^ (^1 - ^ j [n], (2.29) 

are given by 

Ae^iq, t, z) - Aeo(g, t,z)^^{q- l)K'z\-- (^1 " ^) (l ^ ^) K (2-30) 

On taking the results for Qi, Sj, rj, and into account, the Hamiltonians (2.15) of the 
SUSYQM hierarchy can be written as 

K = ^(a^P' + kX^) -SX + a i = 0, 1, 2, . . . (2.31) 

where Oj, bi, Ci are constants independent of S and given by equation (12.36). Going back 
to variables with dimensions, we get 

Hi = hjjjhi — h -miUJ^x^ — qSx + Cihw rrii — — cui — J Uibi u. (2.32) 

2mj 2 Oj ^ 

Note that contrary to the mass and the frequency, the electric field is the same for all the 
partners. 

2.2 Some special cases 
2.2.1 Limit ct 

For small a values, we obtain 

l-2ny^ + 0(a) ^ (^1 - 2s (l - ny/ap^ + 0{a). (2.33) 

Inserting such results in (2.26) transforms the latter into 

I' 

2s 

In the limit a ^ 0, we therefore get 

1 



Ae„~-— + 0(a) where s~l + 0(a). (2.34) 



Aen(S) = -^E^ (2.35) 

which is independent of n and of (3 and is actually the same as in conventional quantum 
mechanics. Wc conclude that the n-dcpcndcnt corrections to an exponential spectrum 
obtained in the general case reduce to n-independent corrections to a quadratic spectrum 
in the special case a = 0. Such a result would also follow from considering equations (2.7) 
and (2.8) for a = 0. 
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2.2.2 Limit /3 ^ 

For small /3 values, we obtain 

g-*^ ~ 1 + 2ny^ + 0{(3) ^ ^1 - ^ j ~ 2(s + nga) + O 
so that equation (2.26) now becomes 



(2.36) 



Ae„ ~ -— ^ + 0(JP] (2.37) 



2(s + nga) 



where 



gc^l + 0{f3) sc^^a + Jl + ^a^ + 0{P). (2.38) 



In the limit /5 — > 0, we therefore get corrections 

-2 



Ae„(a,£) = 



(2.39) 



to a quadratic spectrum e„(«). The latter can be obtained from equation (13.6) by substi- 
tuting a for (3. We conclude that the essential feature of the general case, namely the n 
dependence of the correction terms, is already present when there is only a nonzero minimal 
uncertainty in momentum. 

The result (2.39) may also be derived directly from SUSYQM and shape invariance 
without resorting to a limiting procedure. Going back to the factorization and shape 
invariance conditions given in section 2.1 and setting (3 — therein, we are only left with two 
parameters and r^, satisfying the conditions Sj+i(si+i — a) = Si{si + a), rj+i = rjSj/sj+i, 
since gi = g = 1. This leads to Si = s + ia and n = —E/si, from which equation (2.39) can 
be easily obtained. The fact that r.^ is not independent of i, contrary to what happens for 
a = 0, is clearly responsible for the n dependence of Aen{o.,£) shown in (2.39). 

2.2.3 Case a = /3 ^ 

In I, we showed that whenever the two dimensionless deforming parameters a and /? are 
equal (which means that there is a specific relation a — m^uj'^f3 between the original pa- 
rameters a, (3 of equation (2.3)), the harmonic oscillator Hamiltonian |(-P^ + X'^) reduces 
to a gr-deformed harmonic oscillator one 

/iosc = i(g + !){&,&+} (2.40) 
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where 

satisfy the relation 

bb+ - qb+b = I (2.42) 

with q = (1 + a)/(l — a). 

In terms of the operators (2.40) and (2.41), the Hamiltonian (2.4) can be rewritten as 



^ = ^osc-2V^ + l^(^ +^)- (2-43) 

Since for a — P, we get ^ — 1, t — and u — 2/-\/l — a — \j2{q + 1), it follows from 
equation (2.26) or (2.27) that the correction term to the eigenvalues (13.11) of /iosc reads 

Ae^{q,S) = —^q-^ (2.44) 
q + l 

or 

Aen{q,z)^-^K'z'q-- K = ^{q + 1) z^2£{q + ir'/' (2.45) 
which is again n dependent. 

2.3 Hamiltonian eigenvectors 

Before going to the general case 7^ a 7^ /3 7^ 0, it is worth considering the special case 
a — P ^ 0, ior which the Hamiltonian eigenvalues and eigenvectors depend only on two 
parameters instead of three. 

2.3.1 Case a = f3 ^ 

In terms of the g-deformed boson creation and annihilation operators b^, b, defined in 
equation (2.41), the operators B'^^gj, Si,ri) corresponding to the Hamiltonian (2.15) of the 
SUSYQM hierarchy, assume the simple form 

B+{q, Zi) = ^Kq^/'{b+ - z,) B-{q, z^) = - z,) (2.46) 

with 

z, = zq-\ (2.47) 
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In the (g'-deformed) Bargmann representation of the operators 6^, b, associated with the 
corresponding g-deformed coherent states [26, 27], they are represented by 

Zi) = ^Kq'^'i^ - Zi) B'iq, z,) = ^Kq'/'iV, - z^) (2.48) 

where ^ is a complex variable and Vg is the g-differential operator defined by Vgipl^) — 

[m)-m]/[{Q-m]- 

Here it should be stressed that although the electric field remains constant in the Hamil- 
tonians hi of the SUSYQM hierarchy, the corresponding parameter Zi in the first-order 
operators B^{q, Zi), involved in the factorization of hi, is i dependent, as shown in (2.47). 

The ground state of h is represented by a function ipolq, z;^), which is a normalized 
solution of the equation 

B-{q, z)Mq, -2; = or Vg^Q, = ^Mq, 0- (2-49) 

It is therefore the Bargmann representation of a g-deformed CS \z)q with real z and it is 
given by [27] 

MQ,^;0-J^o{q,z)Eg{zO K{q,z) = [Eg{\z\')]-'/' (2.50) 

where Eg{^) — Z^^oC^/Ng' is a g-exponential. 

The normalized excited state Bargmann wave functions can be determined recursively 
through the equations 

'^n+i{q, z; = [en+i(q, z) - eo(q, z)]~^/^B^(q, z)jjjn{q, Zi, 

- + "(e--)^n(^'J^) (2-51) 

where n = 0, 1, 2, It is straightforward to show that they are given by 

Uq: ^; = K{q: z) n U - ^] Eg (-^A (2.52) 



where 



K(q,z) = {[n],!(g-^"+^(l - q)z';q)J~"^"K (^q, ^) (2.53) 

and n = 1, 2, . . . . In (2.53), the symbol (a; g)„ is defined as in equation (IA.5). 

For z ^ 0, the wave functions ilJn{q,z;^) reduce to the Bargmann representation 



Vn{q',0 = ^^/yb^W- of the n-g-boson states \n)g = {b'^)"'/ y[n]q\ \0)q, which are the eigen- 
vectors of hose- 
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2.3.2 General case / a 7^ /3 7^ 

In the general case of unequal parameters a and g-boson creation and annihilation 
operators satisfying equation (2.42) are defined by 

The operators B^{gi, Si,ri), corresponding to hi, can be written in terms of them as 

B+iq, U, Zi) = -^Kq'/^b^ - t.b - z,) B-{q, U, z,) = -^Kq'^'ib - tM - z,) (2.55) 

and are represented by 
-B+(g, Zi) = -^Kq'/\^ - UV, - zi) B-{q, z^) = -^Kq'/\V, - - Zi) (2.56) 

in Bargmann representation. In equations (2.55) and (2.56), Zi is given by 

zq-^{l-t)(l-^] ' (2.57) 



and reduces to the value given in (2.47) ior a — P or t — 0. 

The ground state Bargmann wave function now satisfies the first-order difference equa- 
tion 

V^Mq, t, ^; = {tC + z)A{q, t, z; 0- (2.58) 

Prom the latter, it is clear that for large values of |^|, the behaviour of '^^{q^t, z;^) is the 
same as that of ^o(?, ^; 0> corresponding to £^ = 0. It therefore follows from I that equation 
(2.58) has a normalizable solution on the whole complex plane. By using properties (IA.2) 
and (IA.8) of the g-exponential and the g-differential operator, respectively, the latter can 
be written as 

ilJo{q, t, z- = M^{q, t, z)E,{Xi)E,{i^^) (2.59) 

where 

X^-z + A n^lz-A A = J]z^ — (2.60) 

2 ^2 V4g-1 

and JVoi^qjtjZ) is some normalization coefficient. 

On expanding the two g-exponentials in (2.59), ipo{q,t, z;^) becomes 

00 

V'o(g, t, z- i) = Mo{q, t, z) Cn{q, t, z)(pn{q; (2-61) 

n=0 
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where 



Cn{q,t,z) 



n 
m 



^/\n\^. p=0 



n/2] 



n-l 



E hi. E 

ri=2p— 1 r2=2p— 3 

p= 1,2,..., [n/2] 

(oo \ -1/2 

E i> ^) 
n=0 / 



ri-i-2 

• E N. 

ri=2p-2i+l 



rp-i-2 

■ E 



ri,=l 



n 
m 



(2.62) 
(2.63) 

(2.64) 
(2.65) 



-I <? 



Here [n/2] denotes the largest integer contained in n/2 and 

is a g-binomial coefficient. The proof of equation (2.64) is by induction over p. 

Note that A, defined in (2.60), may be real or imaginary according to the relative values 
of the parameters. For instance, for a = P, we get t — and A = z/2 so that X — z and 
fjL — O are real. In such a case, equation (2.59) reduces to equation (2.50). On the contrary, 
for £^ = 0, we get z — and A = \^t/{q — 1) so that A = i^t/ {q — 1) and /i — X* are 
imaginary. Equation (2.59) then becomes 



MQ,t;0=^o{q,t)Eg{u 



t 



t 



■e =N'o{q,t)Eg2 



(2.66) 



where in the last step we used a well-known property of the g-exponential [28, 29]. Equation 
(2.66) coincides with equation (14.12). 

The normalized excited state Bargmann wave functions satisfy the recursion relation 



tpn+i{q,t,z;^) = <[n + l]Jl 



^n+l 



1 + {q-l)z'q 



2-n-l 



t 



in+l 



-1/2 



(2.67) 



(2.68) 



X (^ - Wg - z)ipn{q, ti, zi, n = 0, 1, 2, . . . 

where ti — t/q and Zi — {z/q){l —t){l — t/q)~^. 
The solution of equation (2.67) can be written as 

V'n(g, t, Z; i) = Nn{.q, t, z)Pn{q, t, Z; ^)Eq{XnOEq{llnO 

where A„ and /x„ are given by equation (2.60) with t and z replaced by t„ and Zn, respectively, 
Pn{q, t, z; ^) is an nth-degree polynomial in ^, satisfying the relation 

P„+l(g, t, Z; ^) = (^ - z)Pn{q, h, Zi, - t{tn+l^ + Zn+l)Pn{q: h, Zi, q^) - WqPn{q, ti, Zi, ^ 

(2.69) 
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with Po{q,t, z;^) = 1, and J\fn{q,t,z) is a normalization coefficient fulfilling the recursion 
relation 



J\fn+i{q,t,z) 



l + {q- l)z'q-'^-' 1 - 



t 



-1/2 



xNn{qM.Zi). (2.70) 
It can be shown that for the first few n values, the polynomials Pnl?, ^, z] ^) are given 



by 



Pi{q,t,z;0 
P2{q,t,z;0 



-In 



(2.71) 



1 - 



1 - 



(2.72) 



-1 



P3{q,t,z;0 = 



+ [3L 



([2], + 5)-([2], + l)t 



go 



3n -1 



3 



5. 



go 



(2.73) 



For z — 0, they reduce to the corresponding polynomials obtained in (14.26) - (14.28), while 
for t — 0, they give back the polynomials multiplying the g-exponential in equation (2.52), 
as it should be. 



2.4 Comparison with a previous work 

Prom equation (2.9) and equation (2.55) corresponding to i = 0, it follows that the Hamilto- 
nian h, defined in (2.4), is a Hermitian form bilinear in the g-boson creation and annihilation 
operators b^, b, satisfying equation (2.42) with g > 1. The spectrum of the most general 
abstract Hamiltonian of such a type was investigated some years ago [25]. It is therefore 
interesting to see how our results compare with those previously obtained. 

In [25], it has been shown that there exist two factorization schemes leading to two 
possible branches of the discrete spectrum for the considered Hamiltonian. In our case too. 
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the existence of such schemes has been noticed: the parameter k = g/s indeed satisfies 
some quadratic equation, of which we have only kept the root with a plus sign in front of 
the radical, as shown in (2.11) (see also [14]). Our motivation for eliminating the other 
root with a minus sign in front of the radical has been purely physical. We indeed con- 
sider as previous authors working in the field of very small quadratic corrections to the 
canonical commutation relations leading to nonzero minimal uncertainties in position and 
momentum [4, 5, 6, 7, 8, 9, 10, 11, 12, 13] that our theory should give back the standard 
results when the deforming parameters a, /? go to zero. The other solution for k appears 
to be connected with a so-called classically singular representation (see e.g. [30]), which is 
regarded in such a theory as unphysical. 

Focusing now our attention on the first factorization scheme of [25], we notice that our 
result for the energy spectrum, given in (2.24), (2.25) and (2.27), can be retrieved from 
equation (20) of [25] after identifying the parameters ckq, Pq and 70 of equation (23) in the 
same reference with -fC/v^, —Kt/\/2 and —Kz/\/2, respectively. 

In section 2.3, we have provided a thorough study of the corresponding Bargmann wave 
functions, which in contrast with what is claimed in [25] does not happen to be very simple. 
It indeed turns out that whilst the ground and first-excited states (see (2.59), (2.68) and 
(2.71)) arc correctly given by equations (38) and (39) of that reference after identifying the 
parameters as above-mentioned and z with such is not the case for the higher-excited 
states, which as immediately seen from (2.72) and (2.73), do not factorize as stated in [25]. 
Furthermore, our approach has directly led us to a compact expression for the ground-state 
wave function (2.59) as a product of two g-exponentials instead of the infinite product of 
quadratic factors displayed in [25]. The former can be easily expanded into the latter by 
employing well-known properties of q'-exponentials [28, 29]. 

3 Z>-dimensional harmonic oscillator 

Let us now consider the D-dimensional harmonic oscillator problem described by the Hamil- 
tonian 



where the position and momentum components Xi, Pi, i — 1, 2, . . . , D, satisfy modified 
commutation relations of the type [6, 8, 9] 




(3.1) 
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\Pi,Pj] = 

[Xi,Xj] = Ih l+^p2 [PiXj-pjXi) 



(3.2) 



where > 0. Such commutation relations imply isotropic nonzero minimal uncertain- 
ties in the position coordinates Axoi = Axq = ^^^^9-1-^, but none in the momentum 
coordinates, which are simultaneously diagonalizable. 

In the momentum representation, the operators pi become multiplica- 
tive operators, while the operators realized as differential operators 



in 



[1 + (3p'^)d/dpi + f3 PiPjd/dpj + ^Pi 



Here 7 is an arbitrary constant, which does 
not appear in the commutation relations (3.2) and only affects the weight function in the 
scalar product in momentum space 



{f\9) 



where 



a — 



1 + (/3 + 0)p^ 



(3.3) 



(3.4) 



Since the Hamiltonian (3.1) is rotationally invariant, its momentum space eigenfunctions 

can be expressed as a product of a D-dimensional spherical harmonics Yi^_-^...i2ii{^) and a 
radial wave function Rni{p) (where / = Id-i)- In terms of a dimensionless variable P = pa/h 
and dimensionless parameters /3 = Ph^/a^, P' = p'h^/a^, 7 = ^h^/a^, where a = ^Jh/inuu), 
the radial differential equation [9] can be written as 

2 



D - 1 



+ ((D-l)/? + 27)P 



+ 



1 + P^L' - 7(D/3 + /3' + 7) P2 Lr^^(p) = CnlRnliP) 



(3.5) 



where is the eigenvalue of the square of the D-dimensional angular momentum 

L'^^l(l + D-2) Z = 0,1,2,... 



(3.6) 



n is the radial quantum number and 



/(P) = l + /3oP' Po^P + P' e 



nl 



Enl 

fyjj 



(3.7) 
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The first-order derivative in (3.5) can be efiminated by setting 



The resulting equation reads 



(P) 



where 



eni 



1 

2 



^ (L>-3)(L>-1)= / + 



4 

en; - P 



D-3 



/ + 



i:) - 1 



+ -{D-1W 



(3.8) 
(3.9) 

(3.10) 

(3.11) 
(3.12) 
(3.13) 



Prom (3.8), it follows that in the space spanned by functions x(P), the scalar product (3.3) 
becomes 



ix'lx) = I 

Jo 



~ dP 



f(P) 



x'*{P)x{P)- 



(3.14) 



One should note that the arbitrary constant 7, appearing in the momentum space realization 
of Xi and in the scalar product (3.3), is absent from equations (3.9) - (3.14). Prom the very 
beginning, it is therefore obvious that neither Cni (hence e„/) nor XniiP) can depend on it. 

3.1 Energy spectrum 

Let us first prove that /i(^) can be factorized as 



where 



h^'^^B+{g,s)B-{g,s) + io 



BH9,s) = \Uf{P)4T, + 9P-^ 



(3.15) 



(3.16) 



^y^'dP ■ P^ 
and eo is the factorization energy. Here g and s are assumed to be two positive constants 
that are functions of I and of the parameters /3, (3' of the problem. 
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On inserting (3.16) in (3.15) and comparing the result with (3.10), we get the three 
conditions 



s{s-l) = a^'> (3.17) 
g{g-Po) = b^'^ (3.18) 
eo = gs+^ig + M- (3.19) 



Their solution is given by 



s = l + ho-l) (3.20) 

g = ^/3o + A« A« = ^l + /32L2 + i(D/3 + /302 (3.21) 
60 = 1/3o(2/ + 2^-^) + (/ + ^)a«. (3.22) 

In the limit /3, — > 0, we get s = Z + |(£) — l),g'— >1, eg— + which correspond to 
the usual factorization for the D-dimensional harmonic oscillator in conventional quantum 
mechanics (see, e.g., [31] for the three-dimensional case). 

The next step consists in considering a hierarchy of Hamiltonians 

hf = S+(^,, s,)B-{g,, s,) + ^ e,- i = 0, 1, 2, . . . (3.23) 

where }iq = 0\ go = g, Sq — s and gi, Sj, i = 1, 2, . . . , are some positive constants. On 

imposing a shape invariance condition similar to equation (2.16), we obtain the set of three 
relations 

Si+i(si+i-l) = Si{si + 1) (3.24) 

gi+i{gi+i - (3o) = gi^gi + Po) (3.25) 

Q+i = gi+i {si+i + 2^ - 9i - 2) + 2'^°^^'+^ ^^'^^^ 
The solution of the first two is given by 

s + i gi^ g + Poi- (3.27) 
The third one, together with equation (3.19), leads to the eigenvalues 

en{g: s) = E = 5 (2n + s + + Pos (2n + ^) + 2/?on'. (3.28) 
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Taking equations (3.20) and (3.21) into account, we can rewrite them as 



2n + l + 



^] AC) 



2 J 



+ ^Po(2n + l + ^^+Po{l + ^ 



2n+^]+ 2/3on^ (3.29) 



or 



~em = + y) + {^N' + DN-L' + D- 



(3.30) 



in terms of / and of either the radial quantum number n or the principal quantum number 
N = 2n + I. Finally, from equation (3.13), it follows that the eigenvalues in the radial 
differential equation (3.5) can be expressed as 



em 



N- 



^) A« + ![(/? + /?') (at + ^) % (/? - /?') (l^ + ^) + . (3.31) 



In the limit P, P' — > 0, we obtain that e^v; and e^i go to ejv = iV + y, which is the 
conventional result. Equation (3.31) shows that for nonvanishing /3, /?', the spectrum of 
H, given by E^i = hLoeNi, is quadratic in N with an additional / dependence absent in the 
conventional case. The values obtained for E^i coincide with those given in equation (57) 
of [9]. 

The results obtained for gi, Sj and q allow us to write the Hamiltonians (3.23) of the 
SUSYQM hierarchy as 



h. 



(0 



fiP) 



_d_ 
dP 



+ 



(0 

p2 



(3.32) 



where 



,(0 



l + i + 



D-3 



l + i + 



D - 1 



(0 



(0 



In the limit (3, (3' ^ 0, a'p remains unchanged while h)i' 1 and c\ 



D 



/?o(/ + ^-l + y)+A« 



(0 



(3.33) 

(3.34) 
(3.35) 



JO 



As in conventional quantum mechanics, the supersymmetric partners coincide formally 
with some radial Hamiltonians corresponding to shifted / values, / + i, where i = 1, 2, . . . . 
It should be stressed that the angular part y;^_j.../2/j(n) of the wave functions being left 
unchanged, supersymmetry only concerns here the radial equation. 
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3.2 Radial wave functions 



In this subsection, we plan to determine the exphcit form of the eigenfunctions Xni{P) of 
h^''\ Since (/i^'^)^ = 0^ with respect to scalar product (3.14), such eigenfunctions satisfy 
the orthonormality relation 

I -J^^X*nn{P)Xnl{P) = 5n',n. (3.36) 

Note that with respect to (3.14), we also have {B~^{g, s))^ = B~{g, s). 

The ground state wave function Xoi{P) = Xo{9: P) of /i^'^ is obtained from the condi- 
tion 

B-{g,s)xo{9:S;P) = (3.37) 

and given by 

Xo{g,s;P)=Mo{g,s)PV{P)]-^^+''°^y^'''°^=Mo{g,s)P'^+^f{^^^^ (3.38) 
Here s and g are given by equations (3.20) and (3.21), respectively, A and /i are defined by 

and the normalization coefficient is 

The excited state wave functions Xn,i{P) = Xn{g, s; P), n = 1, 2, . . . , can be obtained 
from the recursion relation 



2 



Xn+ii9,s;P) = [/3o(n + l)(n + A + /x + 2)]-i/2i " - ^ ' V ' 27 * P 

xxn{9i,Si;P) n = 0,1,2,... (3.41) 
where from (3.27) and (3.39), it follows that gi — g+Po, Si — s+1 correspond to Ai = A+1, 

yUl = /i + 1. 

Let us set 

Xnig, S; P) = J\fn{g, s)Pi'''^\z)P'^^'^ [/(P)]4(A+M+1) (3.43) 

where N'ni — ■N'n{g,s) is some normalization coefficient and P2'''^\z) is some (A, /x)- 
dependent, nth-degree polynomial in the variable 
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varying in the range (— 1,+1) (with Pq'^'^\z) = 1). Inserting (3.42) in (3.41) converts the 
latter into the relation 



Pii^\z) = [/3o(n+l)(n + A + /. + 2)]-V2^i(^ 

1 r o^d . 1 p^A+i,M+i)(^). (3.44) 



X - 

2 



-{1 - z"^)-^ + X - H + {X + PL + 2)z 



The differential operator on the right-hand side of (3.44) can be recognized as the 
backward shift operator for Jacobi polynomials, satisfying the property [32] 



[1 - z'^)-^ + X - PL + {X + PL + 2)z 
az 



p^A+i,,+i)(^) = 2(n + l)Piif (^). (3.45) 



Hence the polynomials of equation (3.42) can be identified with Jacobi polynomials, while 
J\fn{g,s) satisfies the recursion relation 

from which we get 

Ki9, s) ^ p(2n + A + , + l)n! r(n + A + ,+ l) \ _ 

' \ r(n + A + i)r(n + + 1) y ^ ^ 

When combining equations (3.42) and (3.47) with (3.8), we finally obtain the same 
result as in equation (58) of [9]. 

3.3 Alternative factorizations 

As a final point, we would like to comment on the factorization of /i^'^ carried out in 
section 3.1. 

Equation (3.23) and the corresponding shape invariance condition show that each el- 
ement of the hierarchy hf\ i — 0, 1, 2, ... , admits two different factorizations: this is 
an example of the so-called two-way factorization [33]. The first-order operators B^{gi, Si) 
involved in the factorization act as shift operators, connecting pairs of cigenstates with the 
same energy, belonging to two consecutive Hamiltonians of the hierarchy. 

The conventional harmonic oscillator radial equation is known to admit a four-way 
factorization [31]: there indeed exists another pair of factorizations associated with another 
hierarchy of Hamiltonians i = 0, 1, 2, As before, the corresponding first-order 
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operators B^\gi, si) connect pairs of eigenstates with the same energy, belonging to two 
consecutive Hamiltonians of the second hierarchy. Since, however, the Hamihonians of the 
two extended hierarchies^ are hnked through the relation /i^^/ = h^P — 2i, the shift operators 
of the second hierarchy also connect pairs of eigenstates with different energies, belonging 
to two consecutive Hamiltonians of the first hierarchy. For such a reason, one can combine 
both pairs of shift operators to construct ladder operators connecting the eigenstates of the 
same Hamiltonian [31] (see also [34]). 

In the case of the deformed commutation relations (3.2), the harmonic oscillator radial 
equation (3.9) also admits a pair of alternative factorizations, where in the counterpart of 
(3.16) we choose g > and s < instead of both g, s > 0. Distinguishing by primes all 
parameters relative to these alternative factorizations from those of section 3.1, we obtain 

s' = -/-l(L>-3) (3.48) 



^/9o + A(^) (3.49) 
-l/?o(2Z + i^-I)-(/ + ^)A(0 (3.50) 



and 



s[ = s' + i (3.51) 
g'i = g' + Poi (3.52) 
~< = 9l{s- + ^-gL,{sU-^ + lpo{4 + sl^ (3.53) 



for i = 1, 2, . . . . For /3, P' 0, we get s' ^ -I - ^{D - 3), g' ^ 1, e'^ ^ -I - ^{D - 4) 
and s- = —l + i — ^{D — 3), g[ — > 1, 2, which agree with the conventional results [31]. 

Proceeding as in section 3.1, we can write the energy eigenvalues ejn{g', s') = J2^o oi 
/i(') as 



2 ) 2' V 2 ) ' \ 2 J \ 2 

+ 2/?om^ (3.54) 

^By extended hierarchies, we mean the sets of Hamiltonians or hf^' obtained by letting i run over Z. 
Some of these Hamiltonians may be unphysical. For instance, in the case of h'^p given in (3.32), unphysical 
Hamiltonians correspond io i < —I. 
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The comparison with equations (3.29) and (3.30) shows that the new quantum number m 
is related to the radial and principal quantum numbers, n and N, through the relation 

m^n + l + hD-2)^^{N + l + D-2). (3.55) 

We have therefore rederived the energy spectrum of /i*^'^ in an alternative way. 

The second Hamiltonian hierarchy /if i — 0, 1, 2, ... , containing (3.10) as its first 
member, is given by an expression similar to (3.32) with 

„f = + + (3.56) 



2 

= 6«+2/3oiA«+/3oV (3.57) 



-.(0' 



/?o(-^ + ^ + l-f)+A« 



(3.58) 



For /3, /3' — > 0, af^' remains unchanged, but bf^' — > 1 and cf^' — > i, so that we recover the 
property h^^- = hf^ — 2i for the Hamiltonians of the extended hierarchies. For nonvanishing 

(3, (3', however, there is no simple relation between them and the only common member 
is h^''^ = Hq^ = Hq^' . Ladder operators cannot therefore be built by combining the two 
different types of shift operators. Such ladder operators have actually been constructed by 
using another method [10] and they appear to be very complicated operators. 

4 Conclusion 

In the present paper, we have continued with the application, initiated in I, of combined 
factorization and shape invariance techniques to quantum mechanical problems in the con- 
text of a theory based on some deformed canonical commutation relations and predicting 
nonzero minimal uncertainties in position and/or momentum. 

To start with, we have determined for the first time the spectrum and the eigenvectors 
of a one-dimensional harmonic oscillator in a uniform electric field S. We have established 
that whenever a ^ 0, i.e., whenever there is a nonzero minimal uncertainty in momentum, 
the correction to the harmonic oscillator eigenvalues due to the electric field depends on 
the energy level and actually decreases when going up in energy. This is true whether there 
is also a nonzero minimal uncertainty in position or no, i.e., whether f3 ^ ot (3 = 0. As 
was shown in 1, in the former case the harmonic oscillator spectrum is exponential, whereas 
in the latter it is quadratic. In contrast, whenever a — and (3^0, the electric field 
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induces a constant shift of the (quadratic) harmonic oscillator spectrum, as is the case in 
conventional quantum mechanics. 

In addition, we have shown that although the electric field remains the same in all 
the supersymmetric partners of our original Hamiltonian, its effect on the eigenvectors is 
rather complicated. The latter have been fully determined ior a — (3 ^ 0, in which case the 
harmonic oscillator with nonzero minimal uncertainties in both position and momentum 
reduces to a ^-deformed harmonic oscillator with g > 1. In the general case 7^ a 7^ /3 7^ 0, 
the ground state and the first few excited states have been given in terms of n-g-boson 
states. In both instances, use has been made of a (g-deformed) Bargmann representation 
of the latter and of g-differential calculus. 

Furthermore, we have shown how our results compare with those of a previous study of 
the spectrum of the most general Hermitian Hamiltonian that is bilinear in g-boson creation 
and annihilation operators with q > 1. 

In the second part of our paper, we have reconsidered the problem of a D-dimensional 
harmonic oscillator when there are isotropic nonzero minimal uncertainties in the position 
coordinates, depending on two parameters /3, P'. We have established that our SUSYQM 
techniques can be extended to deal with the corresponding radial equation in the momentum 
representation. As a result, we have rederived in a very simple way both the spectrum 
and the momentum radial wave functions, previously found through lengthy differential 
equation techniques [9] . Finally, we have commented on various factorizations of the radial 
Hamiltonian and stressed both the resemblances and the differences between the deformed 
case and the conventional one. 

The second part of our paper opens the way to other yet unsolved D-dimensional prob- 
lems. Among these, we may quote that of a D-dimensional harmonic oscillator in a magnetic 
field, which would be in line with the first part of our paper. 
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